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Abstract 

We study the electron dynamics in a 2D waveguide bounded by a 
periodically rippled surface in the presence of the time-periodic electric 
field. The main attention is paid to a possibility of a weak quantum 
diffusion along the coupling resonance, that can be associated with 
the classical Arnol'd diffusion. It was found that quantum diffusion is 
possible only when the perturbation is large enough in order to mix 
many near-separatrix levels. The rate of the quantum diffusion turns 
out to be less than the corresponding classical one, thus indicating the 
influence of quantum coherent effects. Another important effect is the 
dynamical localization of the quantum diffusion, that may be com- 
pared with the famous Anderson localization occurring in ID random 
potentials. Our estimates show that the quantum Arnol'd diffusion 
can be observed in semi-metal rippled channels, for which the scatter- 
ing and decoherence times are larger than the saturation time due to 
the dynamical localization. 



1 Introduction 



In this paper we consider the quantum dynamics of a particle in two-dimensional 
rippled channel subjected to time-dependent electric field. In the absence of 
the electric field, the particle motion can be either regular or weakly /strongly 
chaotic, depending on the model parameters. In the case of regular motion, 
the linear response to an external electromagnetic field can be described by 
the famous Kubo formula. In the opposite case of a strong chaos one can use 
the approach recently developed in Refs. PQI2]- The situation with weak chaos 
is quite peculiar and needs a special approach. Of specific interest is the case 
when chaos is extremely weak, and occurs in narrow regions of slightly de- 
stroyed non-linear resonances, thus leading to the so-called Arnol'd diffusion 
along the resonances j3 Hj. Our further analysis enlightens the connection 
between quantum and classical mechanisms of the Arnol'd diffusion on the 
rippled channel model. The results presented here are based on the theory 
developed in Ref. and expand those preliminary reported in Ref. 0. 

The dynamical chaos in classical Hamiltonian systems is related to the de- 
struction of separatrices of nonlinear resonances [7 a . For a weak interaction, 
chaotic motion occurs only in the vicinity of separatrices of the resonances. 
On the other hand, inside the resonances the motion remains regular, in 
accordance with the Kolmogorov-Arnol'd-Moser (KAM) theory (see, for ex- 
ample, Ref. If the number N of degrees of freedom larger than 2, the 
KAM surfaces do not separate the stochastic layers, therefore, they form a 
stochastic web that can cover whole phase space of a system. Thus, if tra- 
jectory starts inside the stochastic web, it can diffuse throughout the phase 
space. Such a diffusion along stochastic webs was predicted by Arnol'd in 
1964 and since that time it is known as an universal mechanism for insta- 
bility and chaos in generic nonlinear Hamiltonian systems with N > 2 (see, 
for example, review [3] and references therein). 

The chaotic dynamics of a quantum systems under time dependent peri- 
odic perturbation was studied in a series of works. The authors of Ref. pQ 
in the frame of random matrix theory have investigated the spectral prop- 
erties of evolution operator in the generic system with a time-dependent 
Hamiltonian. The statistical properties of quasienergy spectrum, the local- 
ization of the eigenstates of evolution operators as well as the process of 
saturation of the energy absorption in the external periodic field was consid- 
ered. In the paper [2| the response of the quantum dot electron system to a 
periodic perturbation was studied analytically in terms of zero-dimensional 
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time-dependent nonlinear sigma model. The quantum correction to the en- 
ergy absorption rate as a function of the dephasing time was calculated. In 
particular, it was shown that the dynamical localization corrections simi- 
lar to the (i-dimensional weak localization corrections to conductivity if the 
perturbation is a sum of d incommensurate harmonic functions. A typical 
application of these results would be to the response of a different elec- 
tron mesoscopic systems (quantum wells, wires and dots) to electromagnetic 
radiation. Experimentally the response of 2D electron gas in quantum dot 
formed in GaAs/AlGaAs heteroj unction to electromagnetic radiation was in- 
vestigated in [3110], where the effect of absorption saturation in open chaotic 
quantum dots was observed in particular. 

The chaotic nature of the free particle dynamics in a rippled channel 
(without external electromagnetic fields) has been investigated in Refs. fTJ 
H21 E3 HI] 5 both in classical and quantum models. In particular, in Ref . [TT] 
the transport properties were considered in a ballistic regime. The energy 
band structure, eigenfunctions and density of states have been analyzed in 
Refs. [121 EH- The structure of quantum states in the channel with rough 
boundaries, including the phenomena of quantum localization, have been 
studied in Ref. [Hj. In particular, it was found that the eigenstates are very 
different in their localization properties. 

The paper is organized as follows. In Sec. 2 we discuss the classical model 
and the mechanism of the Arnol'd diffusion in the rippled channel. Stationary 
quantum states corresponding to the coupling resonance are studied in Sec. 3. 
In Sec. 4 we build the evolution operator for one period of the external field, 
and discuss its structure. In the same section both classical and quantum 
diffusion coefficients are calculated as a function of the goffer amplitude. 
The nature of the quantum diffusion suppression, as well as the dynamical 
localization, are also discussed here. In last section we formulate main results 
obtained in the paper, and discuss the parameters for which the quantum 
Arnol'd diffusion could be observed experimentally. 

2 Classical Arnol'd diffusion in 2D rippled 
channel 

We study the Arnol'd diffusion in a periodic two-dimensional waveguide, see 
Fig. E It is defined by the upper profile given in dimensionless variables by 
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Figure 1: An example of the electron trajectory in the rippled channel. 

the function y = d + a cos x and by the low profile that is assumed to be flat, 
y — 0. Here d is the average width, a is the ripple amplitude, and the goffer 
period is equal to 2tt. In what follows we put the ratio a/dis small, in order 
to avoid the global chaos that occurs for a/d ~ 1. The collisions of a particle 
with these boundaries are assumed to be elastic. 

The Poincare map for the positions and angles (x n , a n ) of the reflections 
from the upper wall is determined by 

f a n+ i = a n — 2 arctan(a sin x n ), , . 

\ x n+ i = x n + tan (2<i + a(cosx n + cosx n+1 )), 

Here x n is the position corresponding to the n-th bounce at the rippled wall, 
and a n is the angle of the particle trajectory makes with the vertical at 

An example of the trajectory described by this map is shown in Fig.|2f a). 
It can be shown that there are many nonlinear resonances arising in this 
model due to the coupling between two degrees of freedom, with the reso- 
nance condition, rj = T x /T y = u y /uj x . Here T y is the period of a transverse 
oscillation inside the channel, T x is the time of flight of a particle over one pe- 
riod of the waveguide, ui x and u y are the corresponding frequencies, and r\ is 
the rational number. It should be noted that in the neighborhood of the reso- 
nances rj = 1/n (for which n — 0, 1, 2, . . .) it is possible to reduce the map ((TJ) 
to the well known standard map [H] with parameter K n = Aad(l + (nn/d) 2 ). 

The mechanism of the classical Arnol'd diffusion in this system is illus- 
trated in Fig. Here some of the resonance lines for the different values 
of rj shown in Fig. H^a), are presented with the use of the u x , c^-plane. The 
curve of a constant kinetic energy E is also shown here, determined by the 



4 



equation 



) 



2 



2E 



"1 + 



(2) 



7T 



m 



where m is the particle mass that we will set to unity. 

As one can see, there are two different types of resonances. The resonances 
with i] < 1/3 are overlapped and in this region the global chaos regime is 
realized. All other resonances are not overlapped, and surrounded by nar- 
row near-separatrix stochastic regions. These resonances are isolated one 
from another by the KAM-surfaces, and therefore, for a weak perturbation 
the transition between their stochastic layers is forbidden. Such a transition 
could occur in the case of the resonance overlap only. In the absence of an 
external field, the passage of a trajectory along any stochastic layer (this 
direction is shown in Fig. for resonance rj = 1 by two arrows) is also 
impossible because of the energy conservation. The external time-periodic 
field removes the latter restriction, and slow diffusion along stochastic layers 
becomes possible. As a result, the particle remains located on the coupling 
resonance, however, with a proportional change of its momentum compo- 
nents. 

The external electric field is given by the potential V(y, t) = —f y(cos Qit+ 
cosf^i), giving rise to two main resonances, u y = Q\ and u y = f^. In order 
to calculate the diffusion rate, we consider a part of the Arnol'd stochastic 
web created by three resonances, namely, by the coupling resonance u x = to y 
and two driving resonances with frequencies Q,\ and Correspondingly, 
we choose the initial conditions inside the stochastic layer of the coupling 
resonance. To avoid a strong overlap of the resonances, however, to pro- 
vide a weak chaos in the near-separatrix region, we assume that the relation 
a/f = 10~ 3 < 1 is fulfilled. 

3 Stationary states at coupling resonance 

In the following it is convenient to rewrite the problem in the curvilinear co- 
ordinates x l for which the both channel boundaries appear to be flat [15J. As 
a result, the covariant coordinate representation of the Schrodinger equation 
has the form 
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where is the metric tensor and g = det(gij). Here we use the units in 
which the Plank's constant and effective mass are equal to unity. Let the 
new coordinates are given by the relations, 

x = x, y = — — ^ (4) 

1 + e cos x 

where e = a/d. Then the boundary conditions read as ip(x, 0) = ip(x, d) = 
and the metric tensor is 



y I ex sin x l+e 2 x 2 sin 2 x J i W/ 

V 1+ecosx (l+ecosx) 2 / 

with the orthonormality condition, 



If the ripple amplitude a is small compared to the channel width d, one can 
safely keep only the first-order terms in e in the Schrodinger equation (JHJ). 
This strongly simplifies numerical simulations without the loss of generality. 
However, one should note that with an increase of roughness in the scattering 
profile y(x) this approximation may be invalid, due to an influence of the so- 
called gradient scattering, see details in Refs.^1]. As a result, we obtain the 
following Hamiltonian [TH] . 

H- H elJ - 1 ( ° 2 ° 2 
" 2 \dx 2 dy 2 

e (~ d 2 n . d 2 d 1 . d \ 

+ - I zcosx-— - — 2ysmx-^-— ycosx— cosx — smx— I . (7) 



2 \ dy 2 dxdy dy 2 dx 

Here and below we omitted tildes in coordinates x and y. 

Since the Hamiltonian is periodic in the longitudinal coordinate x, the 
eigenstates are Bloch states characterized by the Bloch index k. Therefore, 
the eigenstates ip h (x,y) of the total Hamiltonian H can be presented as 
follows, 

^ k (x,y) = J2 ck nm^nm(k,X,y) (8) 



where 



<^ k ^y) = ^/ n+k)x ^{^j L ) ( 9 ) 
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are the eigenstates of the unperturbed Hamiltonian H . Here we normalize 
wave functions to the length L = 2tt of the period in x-direction. In the 
absence of the perturbation (e = 0) the energy spectrum has the form 

KJk) = \ (V + kf + ^) . (10) 

The Bloch wave vector k has continuous values, in particulary, —1/2 < k < 
1/2 in the first Brillouin zone. 

Now we proceed to solve the system of algebraic equations for the coeffi- 
cients c k nm , 

E{k)C nm E nm {Ji)C nrn + e ^ ^ U nm n i m rC n r m r , (H) 



with — oo < n < oo and m = l,...,oo. The matrix elements in (|TT|) are 
defined as 



Unm^n'm 1 — I (Vffc+n'),™') V)' l P(k+n),m,dxdy - 



2 



7r 2 m 2 



(^n',n+l ± ^n',n— l) &m,rr 



+^ ^ f C 1 + 2 ( fc + «))^,n +1 + (1 - 2(fc + n))5 <n _! ) .(12) 

Following Refs. [5], we analyze the energy spectrum in the vicinity of 
the main coupling resonance r) — 1 determined by the condition oj no = u mo 
with uj no = E no+1 (k) - E no {k) = k + n + 1/2 and u mo = E mo+1 - E mo = 
7r 2 (2mo + l)/2d 2 . In a deep semiclassical region for n ^> 1 and m ^> 1, 
one can write u no ~ n and u; mo ~ ir 2 m /d 2 . It should be noted that the 
similar resonance condition can be satisfied for negative hq as well, when 
—no ~ 7r 2 mo/d 2 , that corresponds to the motion in the opposite direction. 
Since below we assume large values, \n\ ^> 1, one can neglect the tunneling 
from the resonances with n > to those with n < 0. Note also that for the 
values of k at the center of the energy band, k = 0, as well as at the band 
edges, k = ±1/2, there are additional integrals of motion due to anti-unitary 
symmetry. For this reason in what follows we consider generic case of other 
values of k. The properties of the energy spectra and eigenstates for specific 
values k = ±1/2 and k = will be discussed separately. 

In the vicinity of the coupling resonance it is convenient to introduce new 
indexes r = n — n$ and p = r + (m — mo). In this notation instead of the 
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system (|TT|) we obtain, 



E(k)c k rp = lpu; mo + ^ + ^(p~r) 2 + kr 

where energy E(k) counts from the level E® omo (k). 

With equation (|T5|l one can investigate quantum states at the coupling 
resonance for different values of the wave number k. For example, by consid- 
ering the states near some positive value Uq and for —1/2 < k < 1/2 (with 
k ^ 0), we obtain that the energy spectrum consists of series of Mathieu-like 
groups. These groups are separated one from another by the energy u no . The 
structure of energy spectrum in each group is typical for any quantum non- 
linear resonance. Namely, the lowest levels are practically equidistant, and 
the accumulation point corresponds to the classical separatrix. However, in 
contrast with the well known Mathieu spectrum, all the above-separatrix 
states are non-degenerate even for the zero goffer a = 0. 

Besides the above discussed part of the spectrum, another series of groups 
occurs in the same energy region, that corresponds to the similar resonance 
with the negative values of n . The eigenstates corresponding to this reso- 
nance are the waves propagating in the negative direction of the x axis. 

In accordance with the spectrum structure it is convenient to characterize 
the states at the coupling resonance by two indexes: the group number q, and 
level number s characterizing the levels inside the group. Correspondingly, 
the energy of each group can be written in the form, 

£U*)="«o(*)9 + ( 14 ) 

where E^ s (k) is the Mathieu-like spectrum for one group. The indexes q and 
s correspond to fast and slow variables characterizing the motion inside the 
classical coupling resonance. Fig. El illustrates the structure of the spectrum 
for one group with q = 0, in dependence on k. The states corresponding to 
no > shown in the figure, have positive values of the derivative dE/dk. 

Let us discuss now the resonance states for specific values k = and 
k = ±1/2. The corresponding states have no classical counterparts because 
the waves propagating in opposite directions are coupled due to the Bragg 
condition. Indeed, at the center and at the edges of the Brillouin zone, where 
the group velocity vanishes, dE/dk = 0, the Fourier series (JBJ) of the Bloch 
function ip k (x,y) contains the terms with positive and negative n equally. 



'rp 
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As a result, the Bloch function ip k (x,y) is the even (and real) or odd (and 
pure imaginary) standing wave. In contrast, the Bloch functions for other 
values of k have no definite parity and can be characterized by a positive or 
negative momentum. 

As one can see, in Eqs. (jHJ) and © for k — it is more convenient to 
use the functions cosnx and sinnx, instead of the exponents e mx . Then, 
the matrix elements U* mn , m , are not different from the analogous ones (|12j). 
and, therefore, the spectrum structure is similar to that discussed before 
for generic values of k. In a similar way, the functions cos(n ± l/2)x and 
sin(n ± l/2)x should be used in the expansion (jHJ) for k = ±1/2. 

The resonance states for k = and k = ±1/2 are either even or odd, 
therefore, the energy spectrum consists of two series of Mathieu-like groups 
corresponding to a specific parity. For k = the above-separatrix levels 
in each group are practically non- degenerated, in contrast with the case of 
generic values of k, see Fig. |3] 

The wave functions of the resonance quantum states have very rich and 
complex structure. To make it clear we plotted the probability distribution 
of the eigenvectors of the Eq. (|13|) with different q and s in the unperturbed 
basis of the Hamiltonian H , see Fig. HJ The probability distribution 
\ip k (x,y)\ 2 in (x, y)-space for the same states are shown at Fig. First, let 
us discuss the results shown in Fig. HJ Here and below, the ground states 
in each group have the number s = 0, and all other states are reordered 
according to the energy increase, labelling by 1, —1, 2, —2, . . . etc. 

So, in Fig. Ufa) one can see probability distribution \c k p \ 2 , corresponding 
to the lowest level in the group with q = 0. In Figs. EKb,c) the probability 
distributions corresponding to the 10th and 20th resonance levels of the same 
group are depicted. To compare with, we show in Fig. EJd) the distribution 
corresponding to the level s = —18 taken from the near-separatrix region. 
All these states are symmetric with respect to n — uq = and m — itlq = 0. 
It is seen that the degree of derealization of these eigenf unctions in the di- 
agonal direction increases with an increase of the energy. Specifically, the 
separatrix state is much more extended in the unperturbed basis, in com- 
parison with the state at the resonance center. As will be shown below, the 
separatrix resonance states, having a maximal variance in the unperturbed 
basis, provide the Arnol'd diffusion. The eigenfunctions corresponding to 
two near-degenerate levels above the separatrix (with s = ±41), are shown 
in Fig. IHe,f). In contrast, the probability distributions corresponding to 
above-separatrix states are non-symmetric (their maximums are shifted with 
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respect to the origin), and have a relatively small variance. 

Fig. illustrates the probability distribution for same resonance eigen- 
states in the transformed (x, ^-coordinates. Note that for small values 
s/d<l the probability distribution in original coordinates has a practically 
similar shape. It is clear that the structure of these states is quite different 
from the unperturbed ones. A regular pattern in Fig. E[a), plotted for the 
state s = 0, resembles the classical trajectory corresponding to the resonance 
center. This orbit starts at point x = —tt, y = and ends at x = tt, y — 0. 
The reflection point is located at x = 0, y = n where the width of the channel 
has a maximum. In accordance with the uncertainty principle, this distri- 
bution has a nonzero width. The probability distributions in Figs. EJb,c), 
corresponding to the 10th and 20th resonance states, also resemble near- 
resonance classical trajectories, however, have a complex internal structure. 
Each state of this kind corresponds to the group of classical trajectories in 
(x, y)-space, for which initial conditions lay at one closed orbit in the phase 
space (sin a, x). One can see that the number of white longitudinal lines 
(minima) are equal to the number of resonance state minus one. Fig. EJd) 
illustrates the probability distribution for one of the near-separatrix states, 
s = —18. Here the pattern also originated from classical unstable trajectory. 
It should be stressed that for all presented states the probability distributions 
mainly have a non-chaotic character. In Fig.[^e,f) we demonstrate the prob- 
ability distributions for above-separatrix states s = ±41. The distributions 
for these states are uniform, indicating that they are weakly perturbed. In 
conclusion, we would like to point out that the electron density distribution 
shown in Figs. EJresemble the electromagnetic field distribution in quasi-optic 
regimes [TT)] . 

Using the same technique we have calculated also quantum states for 
the resonances rj — 1/2 and rj = 1/3. These resonances are somewhat nar- 
row than rj = 1 (see Fig. W[b)) and the number of under-separatrix levels 
are smaller. As an example Fig. H illustrates the probability distributions 
\i/j(x, y)\ 2 for the resonance centers and near-separatrix states for 77 = 1/2 
((a) and (b)) and 77 = 1/3 ((c) and (d)) correspondingly. It is clear that their 
structures are also reflect the character of correspondent classical trajectories 
at these resonances. 
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4 Evolution operator 



In this section we consider the dynamics of a charged particle in the rippled 
channel in the presence of the time-dependent electric field described by the 
potential V(y,t) = —foy(cosQit + cosf^)- We assume that the frequencies 
Qi and are chosen to fulfill the condition u no = (Q± + Q2) /2 in order 
to provide equal driving forces for a particle inside the stochastic layer of 
the separatrix at the main coupling resonance. Specifically, we take, u no = 
400, Q\ = 350, f2 2 = 450, therefore, the period T of the perturbation is, 
T = 7- 27r/fi! = 9 ■ 27r/tt 2 « 0.126. 

Since the total Hamiltonian is periodic in time, one can write the solution 
of the non-stationary Schrodinger equation as iPq(x, y, t) = exp (—iegt) Uq(x, y, 
where UQ(x,y,t + T) = UQ(x,y,t) is the quasienergy (QE) function and €q 
is the quasienergy. It can be shown that the QE functions are the eigen- 
functions of the evolution operator U (T) of the system for one period of the 
perturbation. The procedure to determine this operator was described in 
details in Ref. The matrix elements U qtS ^ q i jS i(T) of the evolution operator 
can be calculated by means of the numerical solution of the non-stationary 
Schrodinger equation for different initial states. The eigenvectors A® (in the 
representation of stationary problem with eigenfunctions (jSJ)) and eigenval- 
ues exp (— ieqT) can be found numerically by a direct diagonalization of the 
evolution operator matrix. After that, the evolution matrix U q<Sjq r tS '(NT) for 
iV periods can be written as 



Therefore, the evolution of any initial state can be computed by making use 
of the evolution matrix, 



In order to illustrate the structure of matrix U qtS}q ' }S i(T) we plot the ma- 
trix elements modules in Fig. [7| As one can see, this matrix has global block 
structure, with a cross-like structure at the centers of blocks. When the evo- 
lution operator acts on the initial state C q ^ s = 5 ? , go 5 S)So , the resultant state 
will coincide with the (go, so)-column of the matrix. The edges of "crosses" at 
the block's centers where matrix elements are relatively large, corresponds 



U q ^ANT) = ^^ s ^;,exp(-^ Q iVT). 



(15) 



Q 




(16) 
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to a transition between separatrix states. As a result, the transition be- 
tween such states of neighbor groups (along the coupling resonance) is much 
stronger than those between other states. Note that the separatrix states are 
responsible for the quantum Arnol'd diffusion. 

Our goal is to analyze the dynamics of a particle, initially placed in- 
side the separatrix under the condition that the coupling and two driving 
resonances do not overlap. In Fig. |H] typical dependencies of the variance 

A q = (AiJ) 2 /u;^ o of the normalized energy are shown versus the time mea- 
sured in the number N of periods of the external perturbation, for different 
initial conditions. Here the quantity A q is defined as the variance of a wave 
packet in the g-space, A q = E 9 (<?-<?) 2 E s |C<?,s| 2 > where q = Y.qQ^Es \ G q,s\ 2 - 

The data presented at Fig. |S] clearly demonstrate a different character 
of the evolution of the system in dependence on initial state. For the state 
taken from the center of the coupling resonance, as well as from above the 
separatrix, the variance oscillates in time, in contrast with the state taken 
from inside the separatrix. In the latter case, after a short time the variance of 
the energy increases linearly in time, thus manifesting a diffusion-like spread 
of the wave packet. 

In order to characterize the speed of the diffusion, we have calculated 
classical diffusion coefficients D\, E x i2 and -D1/3 for the resonances 77 = 1, 
r] = 1/2 and 77 = 1/3, correspondingly: 



To suppress large fluctuations of the energy in time and to reveal a stochastic 
character of motion, the averaging in Eq. (|17j) was performed in two stages. 
Here Ej, is the average value of the particle energy over the time Atj cor- 
responding to 10 3 collisions with the rippled wall. The second average has 
been done in the following way. Having the mean value Ei in each interval 
Ati, the difference (E i+ i — Ei) between adjacent intervals was computed, and 
after that the expression in brackets in Eq. (JT7j) was averaged over all these 
differences. 

The quantum diffusion coefficient Di q was calculated only for the reso- 
nance 7] = 1, see Fig. El It was found that the quantum Arnol'd diffusion 
roughly corresponds to the classical one. However, the data clearly indi- 
cate that the quantum diffusion is systematically weaker than the classical 
Arnol'd diffusion. In all cases the ratio between the frequencies Qi and f2 2 




(17) 
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and u mo of the external field was the same, Q\ : u mo : fi 2 = 7 : 8 : 9. 

One should stress that the quantum Arnol'd diffusion takes place only 
in the case when the number M s of energy stationary states in the sep- 
aratrix layer is relatively large. For the first time, this point was noted by 
Shuryak who studied the quantum-classical correspondence for nonlinear 
resonances. In this connection we have estimated the number of the energy 
states that occupy the separatrix layer. We have found that for a = 0.01 the 
number M s of stationary states inside the separatrix chaotic layer is more 
than 10, therefore, one can speak about a kind of quantum chaos in this re- 
gion. On the other hand, with a decrease of the goffer amplitude, the number 
M s decreases and for 1 / sfa ~ 20 it is of the order of one. For this reason 
the last right point in Fig. El corresponds to the situation when the chaotic 
motion along the coupling resonance is completely suppressed by coherent 
quantum effects. This effect is known as the "Shuryak border" [17]) that has 
to be taken into account when considering the conditions for the onset of 
quantum chaos. 

Since the diffusive motion along the coupling resonance is effectively one- 
dimensional, one can naturally expect an Anderson-like localization. Indeed, 
the variance of the QE eigenstates of the evolution operator is finite in the 
g-space. This means that eigenstates are localized, and the wave packet dy- 
namics along the separatrix layer has to reveal the saturation of the diffusion. 
More specifically, one expects that the linear increase of the variance of the 
energy ceases after some characteristic time. 

In order to observe the dynamical localization in our model (along the 
coupling resonance and inside the separatrix layer) , one needs to study long- 
time dynamics of wave packets. Our numerical analysis for large times (see 
curves (c) on Fig. IHJ) has confirmed that after some time t ~ 200T, the 
diffusion-like evolution stops for different values of the amplitude a. For 
larger times, the variance A q starts to oscillate around its mean value. 

This effect, known as the dynamical localization, has been discovered 
in jTSl EJ] for the kicked rotor, and was studied later in different physical 
models (see, for example, Ref. 0] and references therein). One should note 
that the dynamical localization is, in principle, different from the Ander- 
son localization, since the latter occurs for models with random potentials. 
In contrast, the dynamical localization happens in dynamical (without any 
randomness) systems, and is due to the interplay between (week) classical 
diffusion and (strong) quantum effects. 
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5 Summary 



Let us now summarize our main results. First, we have studied stationary 
quantum states for a particle moving inside a 2D rippled billiard. Main 
attention was paid to the states corresponding to the coupling resonance i] = 
1. We have numerically analyzed the structure of these states, by plotting 
their density distributions |c^ p | 2 in the unperturbed basis determined by 
the model with flat boundaries. We have found also that the patterns in 
(x, ?/)-space of the states corresponding to the inside-resonance region can 
be associated with the classical orbits. Also, it was shown that the states 
corresponding to the separatrix region, have a global structure similar to 
the classical orbits as well. It was found that the width of such states in 
the unperturbed basis is much larger than of those belonging to the inside- 
resonance region. 

In the presence of the external two-frequency time- dependent electric 
field, the main interest is related to the dynamics of wave packets along the 
narrow stochastic regions of the corresponding classical system. Extensive 
numerical simulation demonstrate a kind of weak diffusion in the quantum 
model, that can be associated with the classical Arnol'd diffusion along the 
coupling resonance. We have found that the dependence of the diffusion 
coefficient on model parameters roughly follows the classical dependence. 
However, the quantum diffusion is systematically slower than the classical 
one. This fact manifests the influence of quantum effects. 

It should be stressed that the quantum Arnol'd diffusion occurs in a 
deep semiclassical region, specifically for the case when the number M s of 
chaotic eigenstates inside the classical stochastic layer is sufficiently large 
(of the order of 10 or larger). With a decrease of the amplitude of ripple, 
the diffusion coefficient strongly decreases, and for M s < 1 the diffusion 
disappears. Therefore, one can see how quantum effects destroy the diffusive 
dynamics of the wave packets. 

Another manifestation of quantum effects is the dynamical localization 
that persists even for large M s . Specifically, we have observed that the quan- 
tum diffusion occurs only for finite times. On a larger time scale the diffusion 
ceases and after some characteristic time it terminates. This effect is similar 
to that discovered in the kicked rotor model |18j . and found later in other 
physical systems (see, for example, Ref. j3] and references therein). In our 
case the dynamical localization arises for a weak chaos that occurs inside the 
separatrix layer, in contrast to previous models with a strong (global) chaos 
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in the classical description. 

The effects discussed in this communication can be observed experimen- 
tally, e.g. in semi-metal structures. Apart from the experiments with a 2D 
electron gas in GaAs/AlGaAs hetero junctions [HEIL one can consider the 
semi-metal rippled channel with a large number of discrete quantum levels 
in a deep semiclassical regime. For the observation of the dynamical lo- 
calization in a response to electromagnetic radiation the channel, one can 
take d — 1 /im for the width, and I = 2 fim for the goffer period. Then, 
the dimensionless value a = 0.01 corresponds to the goffer amplitude of the 
order 3.2 nm, i.e. to few monoatomic levels. Correspondingly, for the ef- 
fective electron mass m ~ 0.1m e the level number Uq = 400 has the energy 
E no ~ 0.61 eV, and the resonance frequency for the transitions between near- 
est states is u no /2ir ~ 740 GHz. In these units the period of the external 
field is equal to 1.1 x 10 -11 s, and the dimensionless value of the perturba- 
tion /o — 10 corresponds to an electric field E ~ 0.24 V/cm. It should be 
noted that for above parameters the diffusion saturation time, corresponding 
to approximately 200 periods, is of the order of 2.2 ns. It is clear that the 
electron scattering time has to be much larger. The last condition can be 
realized in semi-metal Bi or Sb structures, where the mean free path is of the 
order of 1 mm and the scattering time is about 10 ns for the temperature 
~ 1 K. 
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Figure 2: The resonance structure for the rippled channel model, (a) The 
Poincare section is shown for d = ir and a = 0.01, demonstrating the struc- 
ture of some of coupling resonances, (b) The positions of some coupling 
resonances and two driving resonances (dashed lines) and the isoenergetic 
curve E = 1.6 • 10 5 on the frequency plane are shown, together with the 
resonances widths (by red color). 
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Figure 3: Structure of the energy spectrum in dependence on the Bloch 
number k, for q = and n = m = 400, d = n and a = 0.003, see text. 
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Figure 4: Probability \c^ p \ 2 in grayscale for different resonant quantum states 
(77 = 1) corresponding to the Mathieu-like group q = (see inset), for pa- 
rameters d = 7r, a = 0.005 and k = 0.01, and for n n = m n = 400. (a) lowest 




Figure 5: Examples of eigenstates probabilities s (x, y)\ 2 in the grayscale 
for the same model parameters and quantum state numbers as in Fig. |3J 
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Figure 6: Examples of eigenstates probabilities \ip(x, y)\ 2 (for resonance cen- 
tre and for near-separatrix state) in grayscale for the same model parameters, 
but for the resonances n — 1/2 fa), fb) and n — 1/3 fc), fd). 
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Figure 8: Time dependence of the variance A q for different initial states 
at the coupling resonance for the group with q = 0: (a) lowest level with 
s = 0, (b) above-separatrix level with s = +45, (c) near-separatrix levels 
with s = -21, +22, -23. Here a = 0.01 and f = 10. 
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Figure 9: Classical, Di, Di/ 2 and -D1/3, versus quantum, Di q , diffusion coef- 
ficients in the dependence on the amplitude a of the rippled profile. 
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